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SECTION I 


INTRODUCTION 

-The purpose of this, report is to: 

(1) Give a tutorial presentation of pseudo- random digital 

sequences, their generation and- properties; and - 

(2) Report the results of a study of filtered pseudo-random sequences, 

, . and their statistical properties'.' * - * ** - = *- - _1 -- -** -r 

-.This study_ re suited from a. need for'speeific design information for * 


a‘- pseudo- random signal generator. The generator, to be used in a telemetry 
communications system. test unit, must generate its -pseudo- random : signals 
by -filtering -a long digital "sequence*. -Desired * signal -properties -include: 
re — - ,(l) Approximately -Gaussian- amplitude— probability -density -function; and 
"Signal spectral "envelope -approx-i-raately Chat of Llie 'f 11 Ler being 


-used -in -the generator. -- r — — -- - ' — — - 

— Filtered maximum- length sequences have been used for this, and similar 

12 3 

applications in the past. ’ 5 The results reported were good for low-pass 

filtered sequences when the ratio of digital clock frequency to filter 

cutoff frequency was between fifteen and twenty. However, for higher 

values of this ratio, a definite skewing of the amplitude density function 

4 . 

has been observed. This skewing effect has been confirmed experimentally 
(See Section VI) for* several digital generator configurations. 

R« P. Gilson has reported that the skewing effect can be averted by 

l 

using a digital sequence composed of the Modulo-2 sum of the outputs of two 
maxinum- length digital sequence generators. His paper included the results 
of experiments with nine and eleven stage generators. Near Gaussian shaped 
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probability density functions were obtained for clock to cutoff ratios of 
up to eighty. 

Experiments supporting this study have confirmed Gilson’s results; 

however, it was found experimentally and by computer simulation, that 

sequence generation by summing two maximum-length sequences did not always 

result in a non- skewed amplitude probability density function. To the contrary, 
* * 

in some cases a single maximum- length sequence gave better results in this 
respect than did the sum of two maximum- length sequences. Therefore, finding 
.a relation between the digital sequence generator feedback configuration and 
k- the sequence statistical structure was important to meeting the objectives of 

"ithe .study, . ..5. ;-f :r ; . 5. - 

Th'is. effort includes an investigation „of the phase distribution of the 
?'spectral components-of various binary sequences, and a direct .analysis of the 
Crrala tion between third moment (skewing .indicator) of subsequences of long 
— bi 

r^yields useful .information. when evaluating the quality of a pseudo-random 
1 -sequence. ; . - - - 




The type sequences to be considered in this section can be described 

as "linear shift- register genera ted sequences". The sequences are linear 
$ 

when' the generator feedback opera tions~3rre res trie ted- to" addition, and not 
multiplication. Such sequences obey the laws of linear mathematics such as 
the law of superposition, .Figure 11-1 is an example of a sequence generator. 
Each box represents a digital storage element, and feedforward and feedback 
are accomplished with the use of "exclusive-or", or "modulo-2 half adder" 
gates. The unit is clock driven, and the output sequence can be taken from 
any stage. The generator in the figure is called a "multiple- re turn generator" 
because of the presence of several feedforward and feedback paths. If feed- 




register, then the unit is called a "simple shift register generator" (SSRG) . 


Figure II-2 illustrates a generator of this type. 
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ft..--- I. Figure IX-2 -Simple -Shift Register Generator . 


These digital generators are capable_of _ producing .repetitive sequences 


of length 


L = 2 - 1 bits. 


.CILJJ 


where n is the number of stages in the register. Depending upon the particular 

-fee'dback— configuratxorr; — the~sequencrr'may cycle in less than L clock periods. 

If this is the case, the generator is called a "non -maximum length sequence 

generator*". - . , _ 

'->• 3 ■ ~ ' -- r_ 

The sequence output of a n on -maximum length generator may depend japon 

the initial loading of the generator. In this case, the generator register 

may pass through a few initial stages’;'' to' which it never 'returns, before 

settling into a non -maximum repetitive cycle.. The initial, .output (before 

getting into the repetitive cycle) is called a "transient output". For 

example, consider the generator shown in figure II-3,„ If the initial .loading 




.is 0111 for the four stages respectively, the output is 111001100110..... , 

•for i .ase rcfererci.’ r.c: c-scr.z i. _r.acz.r-, — . — - 

The 1110 is the transient portion of the output and 0110 is. -the cyclic portion. 

\z~ t i T zz. .carz-r- were zrz curr_: secuLr.t-- .c ;= 

In this case, the cycle is of length four, . 

Tij-v rerr-r c-c:s nrtirzzri.' •: z: sire., sr.zcr register gcn_ r ; — r_ . 
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Figure II-3. Non -Maximum Generator with Transients. 

. Figure II-4 is an example of a non -maximum generator with no transients; 
however, the cycle’s sequence depends on initial loading* 



Figure II-4. Non-Maximum Generator with no Transients. 

Table II-l Lists initial -loading and cyclic -sequence -output for.the generator 

in-figure II-4*_, _ . s~zi,z.zz.„. 

Figure II -5 .is a ’’maximum- length generator"^,. The output sequence (except 
jEor .phase reference) .does .not. depend .upon .initial -loading. „ If . tlje .^f^^ence 
(initial loading) were 1111, the output sequence would be 111101011001000. 

This. report. deala-priraariLy-with^simple shif tregister.generators . ^It._ 

c.5 









Figure II-5. Maximum Length Generator* 

can be shown that every multiple- return generator (that has no transients) 
possesses an equivalent 'simple generator. That is 5 a simple generator ex- 
ists that generates exactly the same sequence. 

‘ Before describing a' formal analysis procedure for simple shift register 
generators s several chara'cteri sties of these generators and their sequences 
will be given. 

No simple shift- register generator which utilizes an odd number of feed 
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back taps can produce a maximum length sequence. Any shift register generator 
producing a maximum length sequence would necessarily include the all 1 state 
of the shift register. If the modulo-2 feedback included an odd number of 
taps, the feedback would be 1, and the next state of the register would be 
the all 1 state. In fact, the generator with an odd number of feedback taps 

would "lock-up" in the all 1 condition, and no maximum length sequence could 

« „ 

be produced . 

If an N-stage simple shift register generator has feedback on stages 
n, k, m, ..., the reverse generator will Have taps on stages n, n-k, n-m, 

..... , For example, the generator in figure II-5 has taps on stages 1 and 4 
and generates the sequence 111101011001000. The reverse generator has taps 
on stages 3 and_4_and generates the sequence 111100010011010. The justification 
for this is as follows: 

Represent the state of the register as 

y ^ 2 ) 4 * * 4 

After one clock period, the state of the register becomes 


V V '•* X n-1 


where 




. bj + _x^ + Xj + x = 0, (II-l) 

andJL, j, k, etc. are the stages included in the feedback. Now consider 


the state 


c 1 


X , X . , ... X, 

n’ = n - 1 1 


in_a generator whose feedback is from stages n-i, n-j, n-k, etc. The 
previous state of the register was 




-2’ *•* X 1‘ C 1 


V > r ' _ . - * 

An 1 / /% \ 


X( i) : content vector aunr.n *. clock ncricr 


z 



where 


1 i 

Comparing (II-l) and (II -2) 


c + x,, + x . + x, + x = 0 
j k . n 


notice that 


is the reverse of 


bj - v 


£ ,, x 9 

n-l n-2 


b l’ *1» 


(II- 2) 


(II -3) 


' • • x. , b. 


n-l 


and the second generator produces a sequence that Is the reverse of that 
produced by the first generator. 

An important characteristic of maximum length sequences that Is used 
extensively in the study of their statistical structure is the "shif t-and -add" 
property. A statement of the shif t-and -add property is as follows: 

If a ML sequence is modulo- 2 added with a shifted version of iLself, the 
result is a shifted version of the original sequence. . 

Another similar characteristic of maximum length sequences is the sampling 
property : 

If the sequence is sampled every k bits where k=2®, g an integer, then the 
resulting sequence is a shifted version of the original sequence. These two 
properties will later be verified using the generator characteristic equation. 

The contents of a shift register has been represented as a vector, 

X ss Xj, Xg, ... . (II- A) 

The operation of the generator can be represented as a matrix operating on 
the content vector. The operating matrix, called the "A" matrix, operates 
on the content vector such that 

AX( j) = X(j+i), (II-5) 

where 

X(j): content vector during j ^ clock period 
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Then, from (11-5) 


X(j+1): content vector during j+1 ^ clock period 


n 


where 


Xi<jn, "£i a ^ <3) 


x^ ( j + l) : - i ^ element of the content vector 
during the (j + l) th clock period. 


• (II-6) 


1 . x^(j): element of the content vector during 

the j ^ clock period. 

As an example, consider the generator^ shown in figure II-6. 



Figure II-6, Five Stage Shift Register Generator, 


The A matrix for this generator is 

o ; : i_ o o 

10 0 0 


-A = 


O'T'O O' 
0 - 0-1 0 
o :o o l 


(II-7) 




' u 


vnicr. is i .1 co ic m i t , Ci x t- v&-uo c c 
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This can be verified by noting 


x 1 (j+l) = * 2 (j) + * 5 (j) 
x 2 (j+l) = XjCj) 
x 3 <j+l) = x 2 (j) 

. x 4 (j+l) = x 3 (j) 

x 5 (j+l) = x 4 (j) (II-8) 

The A raatirx can be interpreted as representing the feed from the j ^ 
stage to the i ^ stage. If 

a ij = 1 

-then the j ^ stage is part -of the feed to the i fli stage, otherwise 


th 


-a. . = 0 . 
ij 


-Notice that' the matrix representation can be used to represent both maximum 
and non -maximum length sequences. 


An extension of the operation -of the A matrix yields the relation 

' - " - X(j+k) = A k x(j). - ' : - 

-Now, if the generator is a maximum length type, then 

-X ( j + L) * A L X(j) = X(j) 


(II- 9) 


(11-10) 

-V = 2 n -l. (11-11) 

-Then - : : 

--- -- i='-ac,=_ r: .(11-12) 

'where I is the identity matrix'. Likewise,- if the sequence were non-maximum 


-and of length P where P<L, then — -- -- r — --- -- — 

cau: ic I. 

r The- characteristic polynomial of the A' matrix is 

I* : Ul" ' 

which is the determinate value of A - £ I, j 


( 11 - 13 ) 


i —oc j. ' « a - 



The charac teri s tic„ equa tion of the A matrix is 

|a - i l| = 0. (11-14) 

A formal relationship exists between the characteristic equation of a 

simple shift register generator and its feedback taps. If stages n, k, ra, 

etc. are part of the feedback, then the characteristic equation is 

+ | n ~ k + 4 n " m +...+1=0. (mod ”2) (11-15) 

A short notation for the characteristic polynomial is 

(n, n-k, n-m, ... 0). (11-16) 

As an example of a characteristic polynomial, consider the SSRG of 

figure II- 6. The characteristic polynomial is 

l 5 + 4 3 + 1. (mod 2) (11-17) 

The short form variation of this polynomial is 

~ — 

(5,3,0). 


A short notation for describing the feedback arrangement for a SSRG 
will now be illustrated. If the feedback is from stages n, k, m, etc., the 
short notation for the feedback is 

[n, k, m, ... O]. (11-18) 
The generator illustrated 'in rfigure II -6 has feedback 


using the above notation. 5 V ' - ~ ~ - ' - --- 

y r The sequence generator characteristic- equation - can be related- to the ' 
statistical properties. oiLcthe: sequence, and therefore is important to this 
study. The relation between "the characteristic equation and ‘the feedback-- 
taps is. .derived as f oil-ows;t . - m smlin; m —z-z.= cer.i-e- 

Fxom Algebra.,. wi.th every -polynomial of the form ... -= cc-.r-r.: 

v- 2 :“c-r" t::- ;.s. . A'\.£ n ^+,C. g 1 ' 1 Hi c+_(P-= (mod -2)- - :(-II-l 9) 

^ 1 2 n 


C — 




■■ ** • 


it: 


vectors 



there exist a matrix (called th'e Companion Matrix) of the form 


C = 


1 

0 

0 

1 


0 

1 

0 

1 


0 

0 

1 

1 


c c c 0 c 

n n-1 n-2 n-3 


0 

0 

0 

1 


( 11 - 20 ) 


If A - §1 = 0 is evaluated, then the resulting polynomial has a C matrix 


given by 

C = A R , ” (11-21) 


or the companion matirx is the rotate of the A matrix. Since the A matrix 
Is related to the feedback taps, the companion matrix is also related to the 
feedback taps. Specifically, if tap k is part of the feedback, then 



and 




c i , . C, — 1 « 
-n,n+l-k -_k - . 


(11-23) 


The- characteristic equation will then have a term of the fora _ . 

, 5 n ' k - ..-<11-20 
K 1 j 1C 

tss £ From- the Cal ey-Hamil ton' Theorem^ a' matrix satisfies its own characteristic 

equation/ The- characteris tic equation can therefore be written- 

A- + A n ' k + A n-m . + ... + I = 0 . (mod 2) (11-25) 

= Using- the~ characteristic equation, the sampling principle stated earlier 

can be justified. If x- is an initial content vector the succeeding content 
2 3 4 

vectors-are Ax, A x, A x, ,A x, . . . If every second output bit, is sampled, 
the output sequence could have been produced by a generator with content 
vectors 


U2 



• ' *) 


.2 A* a6 

x, A x, A x, A x, 

2 

and this is equivalent to having an A matrix of form "A ", Now the character- 
istic equation is 

A n + A n ‘ k + A n ~ m +.,.+1 = 0. (mod 2) (11-26) 

Squaring this equation, 

. - A +A + A +.,.+1 = 0 (mod 2) (11-27) 

where all cross -terms drop out because of the properties of modulo-2 addition. 

2 

Notice (11-27) implies that A satisfies the sequence characteristic equation. 

2 

Therefore, two sequencfe generators with A matrix of form "A n , and "A " 

respectively, produce the same sequence since they each satisfy the same 

characteristic equation. Because sampling every second bit can be represented 

2 

as operating with A matrix of form A , then the sampled sequence must be equiva- 
lent to the original sequence. 

This principle can be extended to sampling every k bits, where k=2 8 , g 
an integer. The sampling principle is not only useful in the study of sequence 
structure, but can aLso serve as the basis for the design of high data rate 
error checking systems^ 

The “shif t-and-add" principle can also be justified with the character- 
istic equation. Consider figure II-7 which depicts two -identical sequence 
generators. Each generator has the same A matrix. Assuming that generator A 
was started g bits earlier than generator B, there is a g bit phase shift 

between the two content vectors. The generator A content vector is 

2 3 

x. Ax, Ax, Ax, ...... 


while the generator B has content vector 

A e x, A e+1 x, A g+2 x, A 8+3 x, 


sssbsvw 


13 
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Figure II-7. Partial Correlation of Two 
















where x is the reference content vector. The content vectors of the sum 


register are 

(A + A 8 + 1 )x, (A 2 +A 8 + 2 )x, (A 3 + A 8 + 3 )x, ....... 

or' rewriting, and using the properties of modulo-2 addition, 

(A + A g+1 )x, A (A + A g+1 )x, A 2 (A + A 8+1 )x, 


This is equivalent to- opera ting the generator with an initial loading 

_ . x = (A + A g+1 )x, _ (11-28) 

and. the sequence of content vectors are 

x , Ax.., .A^ , A 3 x , 

The "A" matrix for the new sequence generator is the same as the original 
A matrix, that satisfied^ the original characteristic equation, only the initial 
■content. is. different.- Therefore, the sum sequence, must be a shifted version of 
Jthe. original sequence. 

It is of interest to investigate the number of maximum length sequences 
available from a given number of register stages. It can be shown that every 
generator with n stages is. capable of producing at least one maximum length 

-*■ - — - r 

sequence depending on the feedback programming. Assuming a maximum length 
sequence from an n stage generator is .sampled every k bits (one sample each 
k bits of the sequence), the new sequence will be maximum length if k is 
.relatively prime , to L = 2 n -l, The Euler's Phi function, 0(K), is defined as 
.the. number of integers, less than K and relatively prime_ to X, therefore, 

0(2 n -l) is the number of sampling rates that are relatively prime to the 
maximum sequence length. The set of sequences obtained by sampling at rates 
U (k* relatively prime ~to"2 n -l) includes all possible maximum length sequences 


of length~~li = 2 n -I. However, the set'includes'repetitions of th"e-same sequence 
because if an "A matrix satisfies ‘a "particular characteristic-equation; -so 
does 'A A ""A ,~~etci"'If the sampling rate k^ is related to the sampling 
rate k 2 by 


, , 2g 

*2 e ’ S an inte Ser 
15 



and both are relatively prime to L (k^ relatively prime to L implies k^ 
relatively prime to L)? then both produce the same maximum length sequence. 

The number of times the same sequence appears in the set with population 
d(2 n cl) must be determined in order to find the number of different maxi- 
mum length sequences available from a given generator. Reference 5 
gives a qualitative proof that the sampling rates fall into groups of 
size n, where n is the number of stages in the generator. The number of 
maximum length sequences available then is 

N = . (11-29) 

n 

As an example, consider the case where n = 5. 2 n -l - 31 and the set 
of sampling rates, 1, 2, 3, ... 30 are all relatively prime to 2 n -l. 

Those sampling rates related by the square relation (Mod 2 n -l) are listed 
in sets below: 

I, 2, A, 8, 16 

■ 3, 6, 12, 24, 17 
5, 10, 20, 9, 18 
.7, 14, 28, 25, 19 

II, 22, 13, 26, 21 
15, 30, 29, 27, 23 

notice the size of each set of equivalent sampling rates is n=5. The number 
of different maximum length sequences available from a five stage generator is 

|SC2 5 -1) 30 

• * r = 6 - 

including reverse sequences. 

Equation (11-29) is exact, however, it becomes tedious to evaluate for 
large values of n, and an upper bound on (II -29) is useful. Such an upper 
bound is devised as follows: 


16 



( 11 - 30 ) 


0(2*-l) = 2 n -2 


2 n 

if n ~1 is prime, and this is the largest possible value of 0(2 -1), 


Then 


0( 2 n -l ) 

n 


2 n -2 2 n 

< £ — £ < £_ 

~ ' n n 


- 2 


n 


2 l0£ 2 n 


- 2 n 'l°S 2 n , 


(11-31) 


An upper bound is then 


D = 2”- log 2 n . 


(11-32) 


For example, if n = 16 


U -= 2 l6r4 = 2 1 '~='4096. 


“•'fiost of- the mathematics necessary for~the analysis of the statistical 
properties of maximum length" sequences “has" been presented# The determination 
of the connections -necessary to- produce- a Triaximum length sequence with an n 
stage generator is not central to the purposes of this report; however, a few 
comments-are in order?~-It can be-proveh" that'no- sequence generator whose 
characteristic equation'is a reducible- polynomial can produce a maximum length 
sequence# Some generators whose characteristic equations are irreducible do 
not produce maximum length sequences; however, all generators that are of the 
maximum length type do have irreducible-characteristic-equations. It has 
already' - been proven that only those generators which utilize an even number of 
feedback taps (odd number of' terms 'in the characteristic equation) produce 
maximum length sequences# Therefore the starting-point in the search for 
maximum length sequence generators* is the set’of irreducible polynomials with an 
odd'humber-of terms. This defines the class of polynomials to be- related to the 
statistical properties of maximum lengtfi sequences of modulo -2 sums of maximum 
length sequences,''- ~- vc ' c - l n: - rx - 
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SECTION III - 


SPECTRAL PROPERTIES 

QF 

DIGITAL SEQUENCES 


A characteristic of white noise is that the spectral components have 
a uniform phase distribution because all the signal spectral phases must 
be independently and uniformly distributed- It is of interest to study 
the phase distribution of the spectral components of a filtered pseudo 
random digital sequence. The pseudo random sequences considered here 
are the maximum length sequences produced by simple linear shift register 
generators, or the modulo-2 sum of several such sequences. 

Ta- at“The”phase distribution of the "sequence spectral components can be 
investigated by calculating the Fourier transform of the entire sequence. 

As a first step, <x Lime reference is established . .A maximum length sequence 
“has _ ' ‘ : ' - ----- 

; ; " L = 2 n -l " ~ ' r ' r ‘ (III-l ) 

bits’ in its sequence. If the true reference Ct^O) is chosen as the center’ 


of the m . bit where 
th 



(III -2) 


then -there are m + % bits during the time interval t<0 and the same number 
in the time interval t^Or If the bit interval is T^, the total sequence 
period is 

T - T b (2 n -1) = T b (2m+l). (Ill -3) 

Figure III-l gives^-a graphic presentation of the time frame and time 
-^reference. During" the bit interval* -the time function is defined as follows: 
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BIT VALUE 

TIME FUNCTION 

1 

-1 - ' 

P(t) 

-P(t) 



Figure III -l. Time Frame and Time Reference. 


The time' representation of the sequence is 


m 


- -- “ - S- -A P(t-gTj 

„ g=-ra g b 

where A is the bit value during the g,, time interval. 

- g; L - th - - „r_ _ 


(nr-A) 


Host storage devices available .for use in shift register generators have 
time functions of form 


=z . ^P(t), = u (t+T b /2) - u (t-T b /2) 


where u(t) is the unit step function. 

The Fourier transform of the sequence is 


; m 


FCx(t))'= J* A P(t-gT ) c' jt0t dt 

g o 


- c0 g=-m 


(III -5) 


m 




-fc- -S Ag J P(t-gTr-) €" J - Wt dt-V 
g~-m -« J 


— — v. . 


—XIII -6) 


Now, changing the variable of integration, 


<0 


F<xCt)> » ’ j 4 J P(t) ,-JMtVb) 

. ° .CO 


g--tn 

- S Ac-« T b” 

g--m 6 


’J 


p(t) e”^dt' 


If 


(III -7) 




A repetitive binary sequence (Period T) can be represented as 


where 


P(t) « £ x a 

O'"- 00 


co — — ~ } 
o T 


(III-8) 


(III-9) 


and 


Sin (a co t /2) 
o b 


r a a co t/ 2 
o 


(III-10) 


Substituting (III-8) into (III -7) and integrating, 


m 


Sin (CO T. /2) 


F(x(t)) = s A «-w b B 6te . aw 

g“-m ° CO T/2 o 


(III-ll) 


Reversing the order of summation. 


F(x(t ) ) = £ SinCC °oV 2 ) 6(co.co) 

0£=_« CO t/ 2 ° o=-m 


m 

£ A C 
g 


-jWgT. 


u 


(III-12) 


which is a set of unit impulse functions at u^OdCO^, a=0, +1, 

Each unit impulse represents a harmonic of the time frunction, x(t), and the 
phase of each harmonic (referenced to t=0, the mid point of the hit in 
the sequence) can be calculated from (III -12), 


m 


0(x) = Tan 


£ Sir /— 1 A (A -A ) 

•1 f£l \2 n -l/ 8 -8 


m /2 tT( 

A + £ cos 


g=l 


?Ol\ 

"- 1 / 


(III-13) 


(A +A ) 
g -g 


The A array is, of course, the digital sequence, 

g 

A computer program was written to calculate the phase of each harmonic. 
The program generates the A array from the binary set I -1,1 1 rather than 

g 4 * 

from {l, oj . 
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As an example, consider the generator shovm in figure III-2. 
loading [-1, -l] the se quen ce output i-s - 1, -l t _l_.___Er.om -(III -13) 

of the 0! th harmonic is Sin(^p) (2) 

6 (a) « Tan \ 7 



With initial 
3 the phase 

(III-14) 


mod 2 \ 



Figure III-2. Two Stage Sequence Generator. 



Figure III -3 depects the sequence spectrum. 



Figure III-3. Spectrum of Output Generator [2, 1, o] . 

-- — s — — - 


Notice,-, .that in; figut~e III -3 the third, sixth, etc. harmonic has zero ampli- 
tude. These are the harmonics -t hat 'fall ah multiplies of the clock frequency. 
Table III-l lists the phase of 'each harmonic. The phases listed for the 






third, sixth, etc. harmonics are meaningless because their amplitudes are zero. 


Table III-l.. Phase of each harmonic from generator 2,1,0 


HARMONIC 

PHASE (RADIANS) 

1 

4-17/3 

2 

277/3 

3 

77 

4 

4tt/ 3 + '/r 

5 

27T/3 + n 

6 

7 7 + 7T 

7 

477/3 

' 8 

277/3 


The additional IT added to the fourth, fifth, and sixth harmonics is 

required because of the negative values of the sine envelope for these 

harmonics. A maximum length sequence (sequence length L) has L harmonics 

before the zero crossing of the sine envelope. 

The computer program calculates the plots, the phase as a function of 

harmonic number, and also calculates and plots the statistical distribution 

of phase. The statistical distribution is calculated based on a truncation 

of the harmonics at This is equivalent to filtering with a perfect 

"brick -wall" filter as illustrated in figure III -4, 

The impulse response of the brick -wall filter is 

h(t) = F -1 (H(w)) 3 (III -15) 

or „ 

■T 

Mt)’« J c jWt , (In-16) 

-T 1 

T 


where 


CO a. 
o max . 

L 


(III -17) 


Evaluating (III -16) 

h(t) » Stn(t V T 

tT^/2 




(m- 18 ) 
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Figure III-4, Sine Envelope and Frequency Response of a "Brick-Wall" Filter* 


Figure III -5 is a block diagram of the generation of the time function 


r(t) where 


t- 

r(t) = J x( r )h (t- Od? 

_<x> “ ' 



Impulse 

Response, h(t) 


r(t ) 


Figure III “5 , Generation of r(t) 


(111-19) 



- ■ The pseudo-randomness, qualities of r(t) can be investigated by a study 
of the pha se distribution of its harmonic components, The functi on. r(t), 
is comprised of the harmonics of the sequence truncated at a » and there- 

luaA 

fore the plot of the distribution of phase provided by the computer program 
is the distribution of phase of r_(t). Section VT gives p lots— of— the .phase 
and phase distribution for several maximum length generators, and maximum 
length sum generators, A maximum length sum generator is .one. in which the 
sequence is generated by modulo -2 summing two or more maximum length sequences. 
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SECTION IV 


FILTERED SEQUENCES 

In Section III an analysis of the spectral properties of low-pass 
filtered maximum length sequences was based upon a filter with frequency 
response as/shown in figure IV-1 . 


FREQUENCY RESPONSE 



Figure IV-1, Frequency Response of Type 1 "Brick -Wall" Filter. 



^This filter has an impulse response given by equation (III-18), 


-r - r:: Sin (tT ) , ^ 

^ (t) = (tT T /2) V 

In Section V an" analysis of the statistical properties of a low-pass 
filtered maximum length sequences is based upon a filter with an impulse 
response of 'the- "Brick'-Wal-l JU -type- as- -shown- irv-f igure -IV -2,.. . 
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The frequency response of the Tvne 2 "RnVV -W»li " filter is 


= j h(t) e^dt = — fj WMT b /2 (IV - 15 


Both of these filters are idealizations, no filter has a frequency response 
as shown by figure IV- 1, and no filter has an impulse response as shown in 
figure IV-2. However, each filter serves as an approximation to a realistic 
filter. As an example, consider a first order filter with transfer function 


HCjbJ 


1 



v, s 

(IV-2) 


Figure IV-3 represents the frequency response of the first order filter. 
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Figure IV- 3, Frequency Response of First Order Filter. 


The impulse response is 


’ h(t) « F ' 1 [l!(]W) ] = J' — 

1+ te 

o 

- 27Tio o < t < “ . 

o 

Figure IV-4 depicts the impulse response, h(t). 



Figure 1V-4, 


Impulse Response of the First Order Lowpass Filter. 
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The forms of the frequency response and impulse response, of course, 
depend on filter type and order. For example, as the order of the filter 
increases the quality of the "Brick -Wall" approximations usually increases. 

The primary justification for using the idealized filters in the analysis is 
the degree of correlation between the results of the analysis and experimental 
data obtained, by filtering. 
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SECTION V 


REPEODUCIBILITY op the 

k /rrn TO DHHP 


STATISTICAL PROPERTIES 


This section is an analysis of sane of the statistical properties of 
maximum length sequences. If x(t) is a random digital sequence, the auto- 
correlation function is 
- - - T/2 

Rxx CO = “ J x(t)x(t+Odt, (V-l) 

-T/2 

where the digital sequence is from the binary set j - 1 , l} • If x(t) is 
-arandoa, then the autocorrelation function will have the form shewn in 
figure V-l . -- - 7 - * ~- 



rs;-'- ""Figure V-l. Autocorrelation of a Random Sequence, X(t). 

If x(t) _is a maximum length sequence (sequence length = L) an appro- 
priate autocorrelation function is 

-j ^T/2- ^ ; - - 

. RxxCO = J x<’t)x(t+ r )dt (V- 2 ) 

“T/2 

where I^L T^, and T^ is a hit period . 7 

New assume the sequence, x, is generated by summing (modulo-2) k 
maximum, length, sequences, ^ - _ . - ~ > 7 . rz 



Modulo 2 


x(t) = 2 x (t) 
i=l 


»v 

(V-3) 


when x(t) is from jl, o| , 


(V-4) 


or 

k 

x(t) = n x . Ct) 
i=l 1 

when x(t) is from | - 1 a l|. 

The length of each constituent sequence is 

L = 2 n i-l 


where is the number of stages in the i^ generator. If none of the sequence 


lengths have common terms 


l = n l . 

i=l 


(V-5) 


The autocorrelation function from (V-2) becomes 


. T/2 

Rxx(T) “ - J jq x.(t) iq x.(t+T)dt (V-6) 

-T/2 i=l 1 i=l 1 


or 


T/2 


1 * * 1 

Hxx(t) = “ J* qj x. (t)x. (t+ T )dt. CV-7) 

-T/2 i=l 1 1 

"When using binary digits from the set -1, 1 , the equivalent of the " shift- 
and-add" property becomes the " shift- and-multiply" property, and if 


T s* L. for i=l , 2, 3, . * *k , 1 , 


X.(t)x.(t+T) = x.(t+T. (T)). 

i i ii 


(V-8) 


Substituting (V-8) into (V-7) 


i T/2 1 

1 k 


Rxx( T ) - ~ J" qj X. (t+T. )dt 
-T/2 i=l 


(V-9) 


where eachT^, is some function of t , Since the sequence lengths have no 
common terms, al) possible phase arrangements of the constituent sequences are 
cycled during the total sequence length T, This characteristic of the product 



sequence can be used to change (V-9) to the form 


L 1 L 2 L 3 

Rxx(t 9*0) = ~ 2 2 2 

L Tj=l 1’ 2 =1 T 3 =l 


V K. 

• z, , n x, (r . > 

T, =1 i=l i x 
k 


h h 

“f 2 K.cr.) 2 x CT, ) 

VI V 1 


■ \ 


T =1 
k 


w- 


Now in each sequence there is one more -1 than +1 and 

■vk 


‘ Rxx (t^O) = 


(- 1 ) 


(V-10) 


(V-ll) 


If =0, there is no shift and 

Rxx (t=0) = 1. 


For the case where r=L i for 3<j<k, then 


or 


(V- 12 ) 


Rxx(t) - ~ 

T/2 

/ 

k 

n x. 

1=1 i 

(t)x. ( t+T ) . 
1 

(V-12a) 

— 

-T 12 

i^j 



and multiply property 

1 

Rxx(t) = ~ 

T/2 

/ 

k 

.n x. 

1=1 i 

(t+T. )dt 
1 


• . 

-T/2 

i*3 

- 


, h 

S’ 


k 


^ i 2 

2 

V* 

• » « • 

.n x.(t.) 


L T i=l 

T =1 

T = 

i i=i ii 


1 

2 

k 


* 

I. 

' = i S' 1 

k 

n 

h 

2 x 

. (T . )1 

(V-12b) 

' L T.=l 
J 

L 1=1 

T.=l 

l 



Rxx(T=nL. ) 

J 

u 

-n k " 1 ' 

4 

(V- 12c) 


Figure V-2 given an example of the autocorrelation function as a function of r 
for even and odd values of k. The scale is exageralcd for the purposes of 





illustration. 


The autocorrelation of a repetitive function is repetitive, as demon- 
strated in figure V-2. Notice the similarity of the autocorrelation function 
of a random sequence (Figure V-l) and the autocorrelation of a product maxi- 
mum length sequence (Figure V-2). 

In Section III the Pseudo randomness qualities of filtered maximum length 
sequences was studied from the approach of the phase distribution of the 
harmonics in the filtered sequence. Figure III-4 is the amplitude response of 
the assumed filter. In this section, a different type filter is assumed. 





Figure V-2. Autocorrelation Function of Sum Sequencer from k Generators. 
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An analysis of the randomness properties can be based upon a calculation 
of the third moment of the filtered sequence. Figure V-3 is the impulse 
response of the assumed filter. 



The frequency response of this filter is 


HCjW) = J h(t) £ _ja>t: dt = Sin(Ct ' MT b /2) ^ -joMT b /2^ (V - 13) 

(Ct/2) 




The time output, rCt) (See figure III -5), is 

t 


f. \ 

rvt; 


f* / / * - 

i xv * ) rue -*' ;q r 


M 

= Z x(t-iT, ) 
i=i b - 


(V-14) 


when x is a digital sequence. Equation (V-14) states that r(t) is equal to the 
sum of the last M digits that occurred before t« A set of M digits taken from 
a long sequence is called a ,,N l-tuple". A study of the statistical properties 
of the filter output can then be based on the weight of the M-tuples where 

M-l 


S - 1 j x . 
X i=0 


(V-15) 


is the weight function. 

• is the number of l's minus the number of -l*s in the M- tuple beginning at 

the m^ b term of jxj when the sequence binary characters are j-1, 1 } . 

• 



In statistical theory the following moments are defined, and each indicate 
certain characteristics of the probability density of the variable x. 


MOMENT 

EQUATION 

INDICATED 

CHARACTERISTIC 

1st 

OD 

. J xp(x)dx- 

V* 

"Mean Value 


■- 00 

- 

2nd 

00 

' - J„x p(x)dx - 

N , 

N 

Variance or . 


Ji - - 

CO 

f x 3 p(x)dx“ 

Central Tendency 

y 

3rd 

Skewning Tendancy 





In the above equations' p(x) is the probability density function. - 

Forming a set-of M-tuple weights, S , with L=2 n -1 members for a maximal 

E 6 . 

length- sequence, S will be defined'as follows.' 

P 1 L ’ 1 P“ : - 

S = £ S-. - -- '(V-16) 

L „ 

1 U-U ~ 


The. first moment is 


L-l L-l M-l 

S 1 = £ S m = -f- 23 2) x 

h m=0 m L m=0 x=0 m+1 


M-l L-l M 

-7— Lt L- x , . = — 
L ra+i L 

- i=0 m=0 


(V-17 ) 


From (V-17) the first moment depends ^on the number of bit periods in thefilter 
impulse response, and on the total sequence length, but not on the particular 
maximum length sequence . - 

The second moment is 

y' ~ nr h. r _ 2 

S 'TT f. S m 

- -„m=0 


. L-l M-l' 2 

— — 2) £ m^i 

" iff=0 i=0 


(V- 18) 
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Now, 


M-l 2 

M-2 

M-l 

o X , , 

Z m + x 

= M + 2 Z 

Z 

i=o 

i=0 

+ 

•rl 

n 

*n 


X ,,x 


and 


s 2 = J_ 

b L 


= M + 


L-l . 

M-2 

m-: 

- 2) M 

+ 2 S 

Z 

m=0 

i=0 

r : 

0 M-2 

M-l 

L-l 

-r s 

z 

z 

-- - v.i-0 

j=m 

m=0 


x _..x_ 


K ,.X , . 

m+x m+j 


From the shift and multiply property of maximum length sequences 


* x _l- X , . = x 

m + x m + j m+£_ _ 

where £ is some function of i-j. _ _ 

Equation (V-20) can be simplified using (V-21), 

1 ^ ^ = M-l 

. S = M + -y- Z Z (-1), 

; \ . i=0 j=i + l 

because over the entire sequence length there is one more -1 digit 
+1 digit » 

Now " 


and 


~ -* - - — M-2- M-l - - ■ “ 

Z Z _ .(-1) = - [(M-l) (M-l) -(M-l) 

: - - - 

~-'[m 2 : 3M+2 ]'- [m-i] 


[h 2 -3M + 2 + 2M-2] 

. 

~H - 1 1 ~ 

' h J ’ 




(V-19) 

(V-20) 
_(V-2l ) 

(V-22) 
than - 

I [m-i] 
(V-23) 


“ M [l 


(V-24) 



The second central moment is by definition 

2 1 1 2 
S * = f S (S -S l ) . 
c L _ m 
m=0 


. (V-25) 


Expanding (V-25) 


2 1 ^ 2 1 ^ ^ k -1 ^ 2 

S Z *:fSS^-2S 1 SS+S (sV 

c m~0 ra m-0 ra m=0 


2 12 
S - (S X )^< 


(V-26) 


From (V-24) and (V-26), the central moment is 


s „ 2 - M [1 - ti] - s 


(V-27) 


The characteristic of the second central moment of the weights of M- tuples 
from the maximum length sequence (filter output) is: 


1) S =0 for M=0 and M=L. 
c 


O 

<- / *-» .*-0 d UICUVXIULUU J.U L 

C ^ 


M = o, 


L 

2 ’ 


s 2 = L±1 = 2 n-2 
c 4 


(V-28) 

(V-29) 


when n is the number of stages in the generator. 


3) The value of does not depend on the particular 
maximum length sequence. 

The third moment of the M-tuple weight set, S^, is 

S 3 - £ S 3 

L _ m 
m=0 


x L-l M-l 

S .£ x 3 
L m+i « 

m=0 i=0 


(V-30) 
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M-l /M-l \ / >1-1 M-l 

2 X J = 2 X (m + 2 2 2 

x=0 \x=0 / N i=0 j=x 


' X , .X . . 
. . , m+i m+ ] 

x+1 J 


= (3M-2) 2 x . . 

1=0 m+1 

■ M-3 M-2 M-l 

+ 3 ! 2 2 2 x , .x , . « 

..... . _ . . , m+j m+k 

x=0 j+x+1 k=j+l 


Substituting (V-31) into (V-30), 


* i L-l M-l 

S = ~~~ 2 [ (3M-2) 2 x , . 

L L . ra+x 

m=0 x=0 


M-3 M-2 M-l 


+ 3! 2 2 2 x ,.x , .x , -i. 

■ 1=0 j=m k-j+1 m+1 ra J m k J 

The first double summation of (V-32) can be simplified as follows: 


L-l M-l 


M-l L-l 

rs v* 

Lt L/ X 


Ul £-j X , . — Lt Lt X , . 

111=0 x=0- m 1 i=0 m=0 m 1 


The quadruple summation can be rewritten 


L-l M-3 M-2 M-l 

2 2 2 2 x ,.x . .x , 

r. •_« ---ii i _ --L, m+x m+i m+k 

m=0 x=0 j=x+l k= j+1 J 


M=3 M-2 M-l L-l 

2 2 2 2 x ..x x . (\ 

i=0 j=x+l k=j+l m=0 J 

Over the range of m in (V-34), from the shift and multiply property of 
maximum length sequences, 

X m+i X m+j X m+£ (\ 

where £ is some function of (i-j). 
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Now 


h y _ f N if 

_ n X m + 4' X m + k j . . 

m =0 ^ (-1 1 


. if £=k 
if £=k 


(V-36) 


Let be the number of ordered 3-tuples, (i, j, k), where 


°<iCj<k<M - 1 


such that 


X , .X . =x ,, . 
ra+i m + j ra + k 


(V -37 ) 


Using the definition of B^, equation (V-32) becomes 


„3 (3M-2) (M) j. 3! 

s = - —L xT 


[0 -,3 


(-1) + B 3 L 


(V-3S) 


Now 


0 


CM) (M-l) (M -2 ) , 
3! 


(V-39) 


and 


s 3 = + 3 HL+1) _ 

L X; S> 


The third central moment is 


3 1 k "i jo 

s J s (S -sV 

c L _ m 

m=0 


fv-Anl 


« ~~ Z (S 3 -3S 2 S L + 3S (S 1 ) 2 - (S 1 ) 3 ) 
m=0 


ra 


o 19 i o i o 

« S' 5 - 3S*S + SCS 1 )-" - (S x )^ 


(V-41) 


or 


S 3 4. B 

c L L 3 


2 3 

M M-l 2M 

+ • 7T • 


(V-42) 


Observe in (V-42) if L is large and M ^ L 


6B„ . 


(V-43) 
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The quantity, B^, can be related to the maximum length sequence 
characteristic equation. If (i, j, k) is one of the tiples such that 
0 < i < j <k < M-l and 


X , .X , .= X • , 
m + i m + j ra + k 

then this is equivalent to having -a generator as shown in figure V-4 


(V-44) 


(k-i) < (M-l) 



Figure V-4, Equivalent Generator of the Class. 

The feedback equation is 

[k-i, k-j, o] 

for the generator shown in figure V-4. 

The sequence, x, must then satisfy the trinomial polynomial of the form 

g(Z) * I + Z^- 1 + Z k " x = 0 (V-45) 


where Z is the feedback matrix (A matrix of Section I). 

The sequence must satisfy the original characteristic polynomial equation, 








Multiplying the content vector by A 8 ^ is equivalent to stepping the generator 
through g^ states, and the sequence must satisfy 

A 8 * (I + C^A 1 ) 

n . 

+ A 82 (1+2 C.A 1 ) 

i=l 1 ‘ ‘ 

~ - ’ + * « . . = g(A) . (V-47) 

Then 

g g 

• ' ' f (A) (A 1 + A 2 + = g (A), (V-48) 

or f(A) is a factor of g(A) « Restating, is the number of trinomials of 
power less than M-l that has a factor f(A), the characteristic polynomial of 
the sequence. _ : r _ — 

From equation (V-43) if is small, the third central moment of the 
weights of M-tuples (filter output) is samll. Reference 6 gives a procedure 
for finding the B 0 term for maximum length sequences for given M-tuples. 


.SOTQDUGIBXLITY OF TF’ 
ORIQ®AL FAGF IS'POOh 
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SECTION VI 


APPLICATIONS AS A GAUSSIAN 
NOISE GENERATOR 

The signal obtained by filtering a maximum length sequence has a probability 
density function that is approximately Gaussian. This can be intuitively 
justified by comparing the maximum length sequence to a random sequence of 
binary digits. If the random sequence digits are from the set j-1, l[ , and 
If during the bit interval the probability of a (-1) is B, and the probability 
of a ( + 1) is 1-B, and each bit is independent of all other bits in the sequence, 
then the probability that there will be m^ (-1) bits during bit intervals is 

p ^-V- ■■■_ p m i (i-p)V m i • (vi -i) 

©1 


If (-1) and (1) are equiprobable, then equation (VI- 1) becomes 





^1 l(m 2 -m^)! 


i 





(VI-2) 


From (VI -2), if a long random binary sequence of equiprobable digits is 
examined, the following properties are evident: 

(a) If there are runs of consecutive (-1) or (1) digits there 
should be 2k^ runs of k 2 ~l consecutive like digits and k^/2 runs 
of k 2 + l consecutive like digits. 

(b) In a long sequence the number of (-1) digits should very nearly 


(c) 


equal the number of (1) digits. 

If two long sequences are compared position for position, the 
number of cases where the digits agree very nearly equals the number 
of cases where the digits disagree. This, of course, leads to the 
correlation property. 
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In Section V it was shown that maximum length sequences satisfy the 
correlation property and property (c) above, -A register generating a 
maximum length sequence goes' through all but one of its possible states, 
therefore, (a) and (b) are also satisfied. The justification for calling 
a maximum length sequence a pseudo -random sequence is based on the fact 
that maximum length sequences do satisfy the randomness properties listed 
above for sequence lengths less than L~2 n Tl for a n stage register. 

Section IV considered filtered sequences, and it was shown that the 
.filtering process is equivalent to forming a signal by weighting and adding 
digits of. the sequence. Likewise; it was. shown that the filtering process is 
equivalent to adding harmonic components of the sequence up to the filter cut- 
off frequency. From the central limit theorem, using either view of the 
filtering process,' the output' signal- density- function tends toward a Gussian 
.density, function if the sequence _is_ random, and the number of digits or 
harmonics added Is large ? This is because, for a random sequence, the 
density function of individual digits are .independent, and the sequence harmonic 
components are indep enden tly distributed in phase. 

.The digits and harmonics components of a -pseudo random maximum length 
^sequence are not actually independent, but in many cases they approximate the 
-random. sequence well enough to approach a Gaussian distribution when filtered. 

The conditions from- the central limit theorem can be used co 9 qualitatively 
test how. well a filtered maximum length. sequence approaches a Gaussian 
.distribution. . . _. „ . _ 

t 

-ft. 

.The. tests are: - .... ** 

1) How random is the distribution of the phase components of the sequence 
harmonics passed by the filter. 



2) If Lhe digits are independent, the weights of M- tuples of the sequence 
should approach a Gaussian distribution and the third -moment of the 
weights should approach zero. Therefore, the value of the third 
moment of the weight function serves as a test for the approach to 
a Gaussian distribution-. 

To illustrate the properties of filtered maximum length sequences, the 
* 

results of experimental tests and computer simulation will now be given. 

Figures VI-1 thru VI-6 give experimental plots of the probability density 
function for the filtered maximum length sequence. The sequence generator is 
a twenty' stage unit”, the filter is- third order, and ratios of the clock to 
filter cutoff frequency, - ' 

‘ ' r ~ r " fc/fcl, (VI- 3) 

are 5, “10, 20, 22;5,' 90, 'and' 200- respectively . From Section IV, the number of 
digits added to form the output is approximately B, Therefore, the statistics 


it. , 
^ 1 1 


Mrti i ; 


From Section V,- the third moment of the M-tuple weights indicates skewing, 
and the third moment is related to the number of trinomials of power less 
than M that have the sequence - characteristic polynomial as a factor. From 
"figures VI-5 and VI-6, the" density function skews, therefore, there must be 
many trinomials of power less- than ‘90 and 200 that have the sequence character- 
istic polynomial as a factor.- -' - — 

* 'As another example", “an eleven stage sequence generator and filtering 
'operation was simulated by computer, lhe maximum length sequence generator 
“characteristic polynomial was' " - 




T 11- 9 

a“+a v +i 


(VI -4) 


Figures VI-7 thru VI-21 give Lhe Lime dependent filter output and density 

c£ the : - trr.'-sf cr. * c £ - . - -,^3. - 

function of the waveform. The figures are for filter impulse response periods 


of 10, 12, 16, 20, 30, 40, 70, and 90 bits. Notice that for impulse response 



periods greater than 10 (M-tuple length greater than 10) there is skewing. 

11 9 

This indicates there exists many trinomials that have A +A +1 as a factor 
for the M-tuple lengths given* 

For the purpose of comparing maximum length sequence generators and sum 
sequence generators, a computer simulation was made of two sequences that were 
formed by adding two maximum length sequences. The first sum generator forms 

t 

its output by summing the outputs of maximum length sequence generators with 
characteristic equations (5, 2, 0) and (6, 1, 0). Figures VI-22 thru VI-35 give 
output waveforms and output density functions for impulse response periods 
of 10, 12, 14, 18, 20, 40, and SO bits respectively. It is seen that the 
density functions for this generator, which uses eleven total stages, do not 
skew nearly as much as for the maximum length generator (11, 9, 0). This 
shows that the characteristic polynomial of the sum sequence is not a factor 
of many trinomials with power less than the M-tuple lengths (impulse response 
peri nd ) . ~ ' 

It cannot be concluded that modulo-2 summed sequences always give non- 
skewed density function when filtered. To illustrate this, a computer simulation 
was run for summing and filtering the maximum length sequences wilh characteristic 
equations (5, 4, 2, 1, 0) and (6, 1, 0). Figures VI-36 thru VI49 give the out- 
put waveforms and output density functions for impulse response periods of 
10, 14, 18, 20, 40, and 80 bits respectively. Notice these density functions 
skew more than those for the first sum-type generators, indicating the 
characteristic polynomials of the generators, are factors of many trinomials of 

power less than the impulse response periods. 

To further investigate the randomness properties of the three sequence 
generators used in the example, a computer calculation of the phase distributions 
of the Tourier transforms of the three sequences was made. Figure VI -50 is a 
plot of the phase distribution of the first 300 harmonics of the (11, 2, 0) 
maximum length sequence. Figures VI-64 and VI -74 a re the phase distribution 




of the first 200 harmonics of the sum sequences ((5, 2, 0) + (6, 1, 0)) and 
(5, 4, 2, 1, 0) + (6, 1, 0) respectively. Figures VI-51 thru VI-63 VI-65 
thru VI-73, and VI-75 thru V I-83 give phase density functions for filtered 
versions of each sequence generator- The number of harmonics included is 
listed on each figure, and the relation between the ratio of clock to filter 

cutoff frequencies, B, and number of harmonics included is ' 

* ^ 

— — harmonics included = r- (VI-4) 

D 

where L is the sequence length, and B = ^c/f^, Notice that the phase density 

functions for the sum generator ((5, 2, 0) + (6, I, 0)) has the most uniform 

distributions, and this -agrees with the previous statements about the relation 
* 

between pseudo -randomness quality and the phase density function of the harmonics 

passed by the filter. ~ ' — „ 

7 The phase distributions_i_llustrated v/ere" calculated J>y. computer from 


eq uation- (I X 1-1 317 
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SECTION VII 


RESULTS AND CONCLUSIONS 


. Section II of this report developed. a method of describipg a sequence 
generator by a characteristic polynomial .-equation,... Section III derived the 
.equation for describing the -phase relations -between harmonic components of 
a. digital sequence, and Sections IV, V, and VI describe methods for investigat- 
ing the randomness properties of filtered sequences. Four methods described 
..and illustrated in the report are: 

1) Experimental method.,- investigating .the. statistics 'of an actual 

. filtered sequence,., ^ ara, ra, 

2) -.Computer- simulation with calculation of the time dependent output 

waveform, and the f il ter. pu.tput pro.bahiI.ity density- func.t-ionr,- • - 

3) Computer calculation. of the phase, distribution of the harmonic 


Am « f-o 




4) Calculation of _the_ thijrd -moment of the filter output by finding 
- i= .. „trinomials r that haye_.the sequence'characteris tic -polynomial as- a 


.All. four methods-h&ve certain rdisadvan tages 5 -however, the 'conclusions "of 
.this section .will be phased on .method .four, pxsscat-- - : a --east," 

Equation Cv - 42 ) relate'd the third central moment of the filter output to 

the number Lrinomials at power less than the impulse response period of the 
filter that has the sequence characteristic polynomial as a factor. Equation 
(V-42) shows that the third central moment can oscillate around, 0 for various 
■filter impulse response periods. This result was confirmed by both experiment- 
ation and computer simulation for both maximum length sequences and for 
% * 

sequences formed by summing two maximum length sequences. 
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An ideal sequence generator for a noxs-e source that must produce approx- 
imately Gaussian waveforms for various ratios of the digital clock frequency 
to filter cutoff frequency (impulse response period in bits or M- tuple length), 
would be one whose characteristic polynomial is a factor of no trinomial of 
power up to the maximum impulse response period anticipated. Two maximum 
length sequences that have characteristic polynomials that are factors of no 
trinomial of order less than 500 are 


(23, 18, 12, 6, 0) 


and 

rl'iT-r/- " - (23ri7, 13, 12, 11, 9, 8, 7, 5, 3, 0), 

Both .these generators have twenty-three stages, and the following stages are 
part of the" feedback': ’ - 


£ , - - (23, 18, 12, 6, 0) - feedback tapes on stages 23, 17, 

. l;~H'ran<i 5,(23, 17, 13, 12, 11, 9, 8, 7, 5, 3, 0) - 

feedback tapes on stages 23, 20, id, id, xo," x4, i^:, 
11, 10, and 6. 

Impulse response periods (digital clock to filter cutoff ratio) will be less 
than 500 for the noise source under consideration. Therefore, the filter 
outputs from a noise source using one at the maximum length generators listed 
above will have non -skewing, approximately Gaussian, probability density 



functions . 
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